Bifurcations of two coupled classical spin 

oscillators 

B.Rumpf and H.Sauermann 
Technische Hochschule Darmstadt 
Institut fiir Festkorperphysik 

January 11, 1994 

Abstract 

Two classical, damped and driven spin oscillators with an isotropic exchange 
interaction are considered. They represent a nontrivial physical system whose 
equations of motion are shown to allow for an analytic treatment of local codi- 
mension 1 and 2 bifurcations. In addition, numerical results are presented 
which exhibit a Feigenbaum route to chaos. 
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Introduction 



Within the last few years several studies concerning the dynamical behaviour of clas- 
sical hamiltonian spin systems have been published. The motion of a single spin being 
exposed to an external driving field as well as an anisotropy field was discussed in [1] . 
Questions of integrability for coupled spin systems under the influence of single site 
as well as exchange anisotropy were thoroughly explored in [2] [3] and also in [4]. 
In many cases dissipative effects cannot be neglected, however. In order to investi- 
gate their influence we study as a paradigmatic example the motion of two coupled, 
damped spin oscillators in a static field B z being driven by a rotating field B in the 
x-y-plane. In agreement with microscopic derivations [5] as well as phenomenological 
treatments [6] the damping is described by adding a Gilbert-term to the Landau- 
Lifshitz- equations of motion. Such a situation may also be interpreted in terms of 
two interacting homogenous modes in a magnetic sample. Despite its physical sim- 
plicity - it is in fact the simplest nontrivial situation if the driving field is circularly 
polarized and there are no anisotropy fields - it will be shown to exhibit a rather 
complex bifurcational behaviour. A similar system which however doesn't possess 
permutation symmetry has been studied in [7] with the main purpose of deriving an 
effective g-factor. 



Bifurcation analysis 

For simplicity and to clarify the notation we start with the trivial case of the motion 
of a single spin 

- S = S x H + TS x (S x H) (1) 

where 

H = B (e x cos ujt + e y sin cut) + B z e z (2) 

and r is the damping parameter. Noting that equation (1) conserves modulus, we 
introduce a rotating coordinate system e! x = e x cosujt + e y smu>t, e! y = —e x sin uot + 
e y cos out, e* z = e z , drop the prims after the transformation and project the sphere 
| S \— 1 into a complex plane. This method is well known in the present context, 
see e.g. [8], [9]. We find the following differential equation for the complex variable 

z ~ i+s z ■ 

z = {i-T)B z z-iujz-{i-T)^-{l-z 2 ) (3) 
This equation yields a stable and an unstable fixed point za/b = — ^3fyj^ =t 

\/ ( (t ^r)Bo ~"> 2 + 1 with the eigenvalues £i,a/b = (i - HC 8 * + b oZa/b) - iw, &,a/b = 
£i,a/b of the corresponding Jacobi matrices. As to the projection technique the reader 
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is refered to the appendix where the most general version in our context is given. 
Our basic Landau-Lifshitz-equation for two driven, isotropically coupled classical 
spins are given by 

— S1/2 = S'l/2 X H e ff 1 / 2 + rSi/2 X S1/2 X H e ffi/2 (4) 

where the effective fields 

Heffi/2 = B (e x cosut + e y smut) + B z e z + JS 2 /i (5) 

now contain the interaction terms JSi/2- These equations are invariant under per- 
mutation of the two spin indices. Transformation into a rotating coordinate system 
and stereographic projection as before yields 

zi = (i - T)S 2 2;i - iuzi - (i - T)^(l - z\) + J(i - T)(zi - z 2 ) \ + Zl * 2 

Z 1 + ^2^2 

Z2 = (i- T)B z z 2 - iujz2 -{i- r)^(l - z\) + J(i - V){z 2 - z^-^^ (6) 

A 1 "T Z\Z\ 

This four dimensional system of equations is invariant under time reversal together 
with a simultaneous change of (B z , B ,u, J, T) to (—B z , —B , —u, — J, T). So, having 
obtained e.g. a solution for positive J (ferromagnetic coupling) one gets a correspond- 
ing solution with opposite stability for negative J (antiferro magnetic coupling) just by 
considering a time reversed situation with opposite directions of the magnetic fields. 
Note that one of the five parameters may be eliminated by introducing dimensionless 
quantities. 

It follows immediately from equations (6), that the spins do not interact, if they are 
parallel to each other. So there are two trivial fixed points 

zi = z 2 = z A (7) 

Zl = Z2 = Zb (8) 

in four dimensional phase space. The eigenvalues of the corresponding Jacobi matrices 
are 

A 1/2 = Re£ ± zlm£ (9) 

A 3/4 = Re£ - 2JT ± z(Im£ + 2J) (10) 

where £ = £4 or £b with £a/b = (i — F)(B Z + B za/b) — iv- Note that £a and £b are 
the eigenvalues of the single spin calculated above. At this stage one proves, that for 
Re£ — 2 JT = and Im£ + 2 J 7^ a Hopf bifurcation occurs. The resulting limit cycle 
is stable (unstable) for J < 0(J > 0). In both cases a center manifold reduction of 
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equations (6) which has been peformed in the settings of ref.s [10], [11] respectively 1 
yields for the Hopf normal form 

z = X 3 z + jz 2 z (11) 



(«-r)J 



where z = Z\—Z2 parametrizes the 2 dimensional center manifold and 7 = — (i+ ZA/B z A/B ) 2 ■ 
As usual the unfolding parameter is the real part of A3. The corresponding bifurcation 
lines Hi and H 2 are plotted in the a>J-plane in figures (1) and (2). 
For B 2 = AJ 2 + T 2 B 2 and uj — (1 + T 2 )B Z the eigenvalues A 3 / 4 vanish altogether. 
Performing again a center manifold reduction the following differential equation is 
obtained 

z = A 3 z + 5z 3 + -yz 2 z (12) 

where 

7 = —j(i — T) and 5 = —J(i — r) r 4i ^f . We have a codimension 2 singularity which is 
unfolded by the 2 dimensional complex parameter A3. Further analysis of this normal 
form shows that it describes a situation in which two saddle node bifurcation lines of 
nontrivial fixed points given by | <5A 3 | — | Im(A 3 7) |= , ReA 3 > (for J < 0) meet 
the Hopf bifurcation line (see point A in figures (1) and (2)). 

We now turn to the nontrivial fixed points for which Sfi X i 7^ Sf ^2. For vanishing 
coupling strength there are only two such fixed points being related to one another 
by permutation of the two spins. We shall see that there may be as many as four 
more fixed points for finite coupling. In order to prove this assertion we develop 
an analytical method which permits not only to solve the fixed point problem of 
equations (6) but also to discuss the pertaining question of stability. We introduce a 
new coordinate system e^,e v ,e^ and project stereographically along the e^-axis into 
its eg — e v - plane. The initially unknown directions of these unit vectors are determined 
unambiguously by requiring that e^ and eg are proportional to Si + S2 and Si — S2 
respectively. As a consequence the fixed points are located symmetrically with respect 
to the origin of the eg — e^-plane at the points ±x of the eg-axis. We prove in the 
appendix that the square of the total magnetisation u = \(Si + S2) 2 — {j^) 2 ^ which 
means that < u < 1, is obtained by solving the following third order equation 

T 2 uj 2 (1 + u) 2 (B 2 - AJ 2 u) = [(1 + r 2 )£? 2 + (B z - u) 2 + T 2 B 2 - 4J 2 (1 + Y 2 )u} 2 u (13) 

Having found a solution of this equation, each of which corresponds to two fixed 
points for the spins being connected by permutation symmetry, its orientation i.e. 
that of eg, e v and e^, may be calculated with the help of the the formulas given in the 
appendix. 



1 Incidentally as far as bifurcation theory is concerned we restrict ourselves by just referring to 
these two books. A useful and more recent overview giving a large though nevertheless selected list 
of references may be found in [12]. 
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In order to investigate possible bifurcations connected with these nontrivial fixed 
points as well as their stability we derive first of all an explicit expression for the 
corresponding (4x4 )-Jacobi matrices. Working in the new coordinate system, 
putting Z1/2 = ±x + Azi/2 and considering equations (21)-(24) of the appendix one 
has the four equations 

1 - x 2 

Ai V2 = ±2(a + bi)xAz 1/2 + (a + U)Az 2 1/2 - 2J(i - r) ^ 2 Azi /2 
±J« - r)[(2* + A Zl - Az 2 ) 1 + (±X + A ^ )iTX + A ^ l} 



1 + (tx + Az 2/1 )(tx + Az 2/1 ) 
1 - x 2 

- 2 TO' (14) 

Now any further calculations would become very tedious if one continued to work 
with B ,B Z and uj as independent parameters. This is firstly due to the fact that 
a, b and x are rather complicated functions of the old parameters and secondly that 
one would have to discuss different Jacobians for different fixed points. Therefore in 
accordance with the concluding remarks of the appendix we consider from now on 
a, b and x as independent bifurcation parameters (leaving J and T unchanged) as far 
as analytical considerations are concerned. However, in what follows we switch back 
to the old parameters for physical reasons when displaying our results graphically. 
Linearisation of equation (14) yields 



(15) 



/ 2ax + JT -2bx + J -JTq J \ 

2bx - J 2ax + JT Jq JT 

-JTq J -2ax + JT 2bx + J 

V Jq JT -2bx - J -2ax + JT J 

for the Jacobian. The corresponding eigenvalue equation reads 

A 4 + K 3 X 3 + K 2 \ 2 + K x \ + K = 
with 

K 3 = -AT J 

K 2 = 8x 2 (b 2 - a 2 ) + J 2 {2 -2q + 5T 2 - T 2 q 2 ) 

K x = 32x 2 abJ + 2JT[8x 2 (a 2 - b 2 ) + (1 + T 2 )(q 2 - I) J 2 ] 

K = I6x 4 (a 2 + b 2 ) 2 + J 2 x 2 [8a 2 (l + q) + Ab 2 (q 2 - 1) - 24a6r + Aa 2 T 2 (q 2 - 1) + 
8b 2 T 2 (q + 1) + 8abT(q 2 - 2q)\ 

(16) 

where q = n+^p — 1, i.e. — 1 < q < 1. Using these formulas it is now possible 
to write down the explicit algebraic expressions for the conditions which have to 
be met for the various local bifurcations to occur. They depend on altogether 5 
parameters a, b, x, J, T and read for saddle-node and Hopf-bifurcations K Q = and 
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K\ — KiK 2 K 3 + K K 2 = 0, KiK 3 > respectively. The corresponding bifurcation 
lines together with the Hopf line of the trivial fixed point are plotted in an u — J-plane 
in figures (1) and (2) by evaluating the algebraic expression with a computer. They 
are symmetrical to the axis J = 0. We find two cusp points C\ and C 2 . We recognize 
point A as the codimension 2 point which was already discussed in connection with 
equation (12). There are six nontrivial fixed points inside the triangle A,Ci,C 2 and 
two outside respectively. 

Next we point out that all of the three basic cases of local codimension 2 bifurcations 
of ref. [10], chapter 7, occur in our system. There is a double zero eigenvalue associated 
with linear part 

(Si) (») 

at Fi and F 2 , where saddle node and Hopf lines meet. One finds a pure imaginary 
pair and a simple zero with linear part 
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(18) 



at the points D in fig. (2). The situation of fig. (2) is obtained from that of fig.(l) by 
decreasing the field B . Thereby Fi approaches A and reaches it at a certain value of 
B and coincides with A later on. Before that happens the two points D\ and D 2 , in 
which the Hopf and saddle-node lines meet tangentially, appear simultaneously, they 
persist thereafter. 

An analytic argument for the existence of these types of bifurcation is the following. 
A double zero eigenvalue exists if and only if K = K x = 0. This is true for a = 
0, Ab 2 x 2 + J 2 (q 2 - 1) = at T = 0. Because 



dK dK 
da db 



7^0, 



(19) 



the implicit function theorem tells us, that there are (for x and J fixed) unique 
solutions a(r) and b(T) in the neighbourhood of that point. For the second bifurcation 
in question the two conditions have to be replaced by K = K 3 K 2 — K x = 0, K 2 > 0, 
the remainder goes through unchanged. In a three dimensional parameter space the 
codimension 2 bifurcations are represented by lines instead of points. They are plotted 
in fig. (3). We draw the attention of the reader to the points V\ and V 2 where Fi and 
F 2 meet line A. At present we desist from discussing them any further. 
We finally state that there are two imaginary pairs of eigenvalues and a corresponding 
four dimensional center manifold at point E in figure (1),(2). 
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Chaotic dynamics 

We do not investigate any further the global bifurcations which usually occur in 
connection with the local codimension 2 bifurcations of the last section. Nevertheless 
it may be expected that they are an important cause for chaotic behaviour in our 
system. What we have done however in this respect is to employ the software-package 
AUTO [13] that can follow unstable limit cycles to get at least some numerical survey. 
We found that the stable limit cycles grown out from the nontrivial fixed points for 
J < by Hopf bifurcations pass through a period doubling scenario which ends up 
in a chaotic attractor (fig. (4)). The Lyapunov spectrum as a function of the coupling 
parameter is plotted in fig. (5). The two Lyapunov exponents which are not shown 
are strictly negative. As is characteristic for such cases of broken symmetry [16] 
there are two homologous Feigenbaum routes related to each other by permuting the 
two spins. For J > there is a corresponding repellor. Increasing | J | further we 
get - after a sequence of periodic windows - a reverse Feigenbaum route. It leads 
after a final pitchfork bifurcation to a stable cycle which then persists forever. It is 
remarkable that the symmetric limit cycle which has been created originally, i.e. for 
small | J | by a Hopf-bifurcation of the trivial fixed point and has changed its stability 
character in a subsequent global bifurcation, is the necessary unstable object in that 
last bifurcation. This was demonstrated independently by observing the continuous 
deformation of the final cycle into the initial one along a completely different path in 
parameter space. 

Summary 

We have proven that the four dimensional equations of motion of two coupled damped 
spin oscillators which are driven by a rotating field lead to a surprising variety of 
dynamical and bifurcational phenomena. One finds all the basic codimension 1 and 
2 bifurcations as well as a route to chaos. As our system represents a model for an 
antiferromagnet with constant sublattice coupling it provides an important physical 
example for a case study of this sort of problems. Our results differ strongly from the 
trivial ones for a single spin oscillator as well as those for two coupled but undamped 
spins [2], [3]. The same is true if we compare them with those derived in [17] for two 
torsionally coupled damped and driven pendula. 
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Appendix 

Projecting the following equation 

_ Si = 5i x (B - uj + JS 2 ) + TSi x (Si x(B + JS 2 )) (20) 
along the arbitrary direction e^ — e^x e v into the — e v plane yields 

ii = (i - r)5 € ^^ - i(i - T)B^-^ + (< _ r)B cZl 

. z\-l 1 + z 2 . 
-IU}£ — lo v — iu c zi 

+J(i-D(z 1 -z 2 ) 1 -±^ (21) 
1 + z 2 z 2 

B^, B v , B^ and cu^, u v , lo^ are the new components of B — Bqc x + B z e z and cu = u>e z 
in the directions e^,e v and respectively. Choosing the axes e^,e v ,e^ as described 
in the text the resulting fixed point equations take the form 

1 - x 2 

= Ax 2 ±Bx + C±J(i-T) - (22) 

1 + x l 

Explicit expressions for the complex quantities A,B,C follow by comparison with (21). 
Adding and subtracting these two equations and dropping the trivial solution x = 
leads to 

= Ax 2 + C (23) 

= (S + 2J(i-r)^) (24) 

Obviously equation (24) requires B to be equal to (i — T) times some real number 
which turns out to be B^. This enforces to be zero, i.e. must be perpendicular 
to u — ue z . Introducing the Eulerian angles $,0 and ^ which connect the two 
coordinate systems we have the result O = | and consequently 

Bj: = B cos^cos$ + B z sin^ (25) 
B v = -B sintycos$ + B z cos^ (26) 
B c = B sin<S> (27) 
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The quantities A=a + ib, C=c + id may now be written down immediately as 



a = i(£„-w„-rB € ) (28) 
b = ±(B € -u> € + rB v ) (29) 
c = i(S„-w„ + rB € ) (30) 

d = -I^-^-rs,) (3i) 

Observing that equations (23) and (24) are 3 real equations for the three unknowns 
$, \l/ and x we eliminate $ and ^ and end up with equation (13) of the text. 
Finally, in the context of the discussion of stability of these fixed points we need the 
following formulas: Using (23) to express c and d in terms of a, b and x we solve the 
four equations (28)-(31) for B^, B v , u^, u v in terms of a, b and x as well. Noting that 
(24) is equivalent to B^ + 2Jj^ = and W( = we end up with 

Bt = -f;(l+x 2 ) u^ = -(f + b)(l + x 2 ) 

B v = -^{1-x 2 ) uj v = (£-a)(l-x 2 ) (32) 
B c = -2J*=£ u c = 

As obviously u = yj + uj 2 + uj 2 , B z = B w+ B ""<< and B = sj '(B + B z ) 2 - B 2 we 
have 3 relations which can be interpreted as effecting a transition from the original 
parameters B , B z , wtoa new set of parameters a, b, x leaving J and T fixed, and vice 
versa. 
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Figure captions 



figure 1 

Bifurcations in the 3-u plane with B = 1.2, B z = 1.0,T = 0.1 

Hi and H 2 are the Hopf bifurcation lines of the trivial fixed point. At E there are four 
imaginary eigenvalues. SN are saddle-node bifurcation lines of nontrivial fixed points, 
inside the triangle A,Ci,C 2 there exist six nontrivial fixed points, outside only two. 
Ci and C 2 are cusp points. H 3 and H4 are Hopf bifurcation lines of the nontrivial 
fixed points. F\ and F 2 have two vanishing eigenvalues. 

figure 2 

Bifurcations in the J-u plane with B = 0.4, B z = 1.0, T = 0.1 

Di and D 2 are bifurcations with two imaginary and one zero eigenvalue, Hopf and 
saddle node bifurcation lines are in tangential contact with each other. 

figure 3 

Bifurcations in the J-u-B parameter space with B z = 1.0, F = 0.1. Codimension 
2 points of fig.(l),(2) appear as lines. Note that D is a closed curve. i*\ and F 2 
terminate at Vi andV2 

figure 4 

Feigenbaum-scenario for B = 0.9, a; = 0.9, B z = 1.0, T = 0.1. The Poincare plane is 
defined by S^i = 0, S y i is plotted over the parameter J. 

figure 5 

Plot of the largest two Lyapunov exponents as functions of the parameter J for 
B = 0.9, io = 0.9, B z = 1.0, T = 0.1. The computation was carried out by a QR 
decomposition based method [14], [15]. 
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